Modeling the electronic behavior of 7~LiV205: a microscopic study. 
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We determine the electronic structure of the one- 
dimensional spin-| Heisenberg compound 7-L1V2O5, which 
has two inequivalent vanadium ions, V(l) and V(2), via 
density-functional calculations. We find a relative V(l)-V(2) 
charge ordering of roughly 70 : 30. We discuss the influence of 
the charge ordering on the electronic structure and the mag- 
netic behavior. We give estimates of the basic hopping matrix 
elements and compare with the most studied a'-NaV^Os. 
PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



Low-dimensional transition metal compounds have 
been intensively studied in the past,«ears. Among those, 
the quarter-filled ladder compound™ 2 ] a'-NaV^Os has be- 
come a model substance for the study of spin-charge and 
orbital coupling. The coupling between spin and orbital 
ordering is a central issue in the somewhat moro-|Com- 
plex colossal magnetoresistance (CMR)-materialsB, the 
advantage of a'-NaV^Os as a model compound is the spa- 
tial separation of the two <i X y-orbitals onto two spatially 
separated V 4 5+ -sites, arranged in ladders. These two 
orbitals undergo a charge order transition 4 ! at T c — 34 K 
with 2V 4 5+ ^V 4+ +V 5+ . Simultaneously .-to the charge 
order and the respective lattice ditortionsB, a spin-gapcl 
opens, indicating substantial spin-charge coupling in this 
compoundj-jThe physics of this transition is being studied 
intensivelyufl. 

A much less studied, though not less intriguing sys- 
tem belonging to the same vanadium oxide family is 
7-LiV2C>5. Susceptibility measurements!], as well as 
NMR experimentst2l, on this compound suggest a one- 
dimensional spin-i Heisenberg-like behavior and there 
is no indication of a phase transition at lower tempera- 
tures. To our knowledge, there is no microscopic study 
of the electronic structure of this material discussing the 
magnetic interactions responsible for such behavior. We 
present here a density-functional analysis (DFT) of this 
system and calculate the possible exchange matrix ele- 
ments via the Lineaif-Mufim-Tin Orbital (LMTO) based 
downfolding methodtil and a tight-binding model. 

7-L1V2O5 offers the possibility, due to its close relation 
to a'- NaV^Os, to study the influence of charge order- 
ing on the electronic structure. In particular this opens 
the question on how far the charge ordering and the cor- 
responing crystallographic distortions alter the magnetic 
interactions. We discuss several possible scenarios com- 
patible with the experimental susceptibility for 7-LiV205 
for the underlying magnetic model (i) A zig-zag chain 
model of V(l)-ions, (ii) A double-chain model of V(l)- 



ions, (iii) An asymmetric quarter-filled ladder model. 

Crystal structure.- 7-L1V2O5 has a layered structure 
of VO5 square pyramids with lithium ions between the 
layers. It crystallizesllj in the orthorhombic centrosym- 
metric space group — Prima and has two crystal- 
lographic inequivalent vanadium sites, V(l) and V(2), 
which form two different zig-zag chains running along 
the y axis. Within the layers, V(l)05 zig-zag chains are 
linked to V(2)Os zig-zag chains by corner sharing via the 
bridging 0(1). The existence of two types of V-sites has 
been also verified by NMR experimentsE-X 

In Fig. |l] we show the crystal structure of 7-LIV2O5 
and a'-NaV205 projected on the xz plane. The angle 
between the basal plane of the (nearly) square V(l)Og/ 
V(2)05 pyramids and the x-axis is about +30°/ —30° 
respectively for 7-LiV"205. The basal plane of the VO5 
pyramids in a'-NaV^Og is, on the other hand, nearly par- 
allel to the x-axis. Note that the VO5 square pyramids 
are oriented along x as down-down-up-up in a'-NaV^Os 
while in 7-L1V205 the orientation is down- up-down- uj}_ 

From the structural analysis it has been proposed^ 
that the oxidations of V(l) and V(2) are, respectively, 
V 4+ and V 5+ . The temperature dependence of the sus- 
ceptibility x(T) follows that of a one-dimensional spin- 
i Heisenberg model with an exchange interaction of 
Jexp = 308 K and a gyromagnetic factor g — 1.80. Note 
that depending on the magnitude of the exchange cou- 
plings, Ji -between V(l)-ions in edge-shared pyramids- 
and J;, -between V(l)-ions in corner-shared pyramids- 
along y the system can be treated as Heisenberg zig-zag 
chains (Ji 3> J&) or as Heisenberg double-linear chains 
(Ji -C Jh), compare with Fig. ||. 

A third possible interpretation of the nature of 7- 
LiV^Os compatible with the experimental susceptibility 
relies on the possibility of a partially charge-ordered sys- 
tem i.e. V(l)-sites somewhat closer to V 4+ -oxidation and 
V(2)-sites closer to V 5+ -oxidation. A picture of an asym- 
metric ladder with one electron per V(l)-0-V(2) rung 
would then describe the system in analogy to cv'-NaV^Os 
where the magnetic interactions among the constituent 
ladders are weakll3. In the following, we will investigate 
these three scenarios. 

Band structure.- We have calculated the energy bands 
(see Fig. ||) of 7-L1V2O5 within DFT by employing 
the full-pcitcntial linearized jaiigmented plane wave code 
WIEN97B and by LMToN based on the Stuttgart 
TBLMTO-47 code. We find complete agreement in be- 
tween these two calculations. 

The overall band picture for 7-LiV~205 (Fig. ||) is simi- 
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lar to that of a'-NaV205EJ. The V-3d states give the pre- 
dominant contribution to the bands at the Fermi level 
and up to f» 4eV above it. The lower valence bands 
are mainly 0-2p states and are separated by a gap of 
« 2.2 eV from the bottom of the V-3d bands. There is 
a nonequivalent contribution from the two types of V- 
sites, V(l) and V(2). The four, lowest-lying 3d bands 
at the Fermi level are half-filledEj and are made up pre- 
dominantly of V(l)-3<i and of V(2)-3d states in the ratio 
p(l)/p(2) « 2 : 1 and 3 : 1 depending on the k- values. 
The next four bands less than 1 eV above the Fermi level 
also exhibit a mixture of V(l)-3d and V(2)-3d character. 

The vanadium bands at the Fermi-level are of d xy sym- 
metry (global symmetry) with a certain admixture with 
the d yz state due to the rotation of the basal-plane of the 
V(1/2)C>5 pyramids with respect to the x-axis (see Fig. 
|l|). The degree of admixture for both vanadium types 
is such that the respective V-3d orbitals point -as in a 1 - 
NaV205D- roughly towards the bridging oxygens. In this 
sense, one can regard the electronic active V(l)-3d and 
V(2)-3d orbitals as d xy orbitals rotated around the y-axis 
by angles <p\ — 35° and ip2 = —28° respectively (see Fig. 
!)• 

The most notable difference between the band- 
structure of 7-L1V2O5 and a'-NaV205 is the band- 
splitting at the X and T points in 7-L1V2O5, which is 
absent in a'-NaV205. This splitting is due to the ex- 
istence of two different V-sites in 7-L1V2O5 as we will 
see in the next paragraph. The fact that this splitting is 
big, close to the overall bandwidth, indicates already that 
the microscopic parameters associated with the V(l) and 
V(2) sites must differ substantially. 

Also note that the splitting of the bands at the Fermi- 
level due to the existence of two xy-planes in the crystal- 
lographic unit-cell of 7-L1V2O5 is small and does not oc- 
cur along the path Z-U-R-T-Z. We will concentrate upon 
the discussion of the in-plane dispersion in what follows. 

Microscopic parameters. - In order to determine the mi- 
croscopic model appropriate for 7-L1V2O5 we have an- 
alyzed the band-structure shown in Fig. [| by a (mini- 
mal) tight-binding model with one orbital per vanadium 
site, which generalizes the tight-binding model appropri- 
ate for a'-NaV^C^ to the case of two different V-sitcs. A 
straightforward fit, e.g. by least-squares, is not possible 
since the lower unoccupied V-bands (roughly in between 
0.5 — 1.0 eV), are strongly hybridized with the O-p or- 
bitals and in order to describe the low-energy physics of 
this system these-bands should also be considered. 

In recent yearsell, a new version of the LMTO method 
has been proposed and implemented which is proved to 
be powerful in providing an effective orbital represen- 
tation of the system by integrating out the higher de- 
grees of freedom using the so-called downfolding tech- 
nique. The usefulness of the method lies in taking into 
account proper renormalization effects. However, the 
Fourier transform of the downfolded Hamiltonian to ex- 
tract the tight-binding parameters, results in long-ranged 
hopping matrix elements. 



We have therefore considered a combination of both 
methods. We use only the short-ranged hopping matrix 
elements (see Fig. |^) provided by the downfolding pro- 
cedure as an input for the tight-binding model. These 
matrix elements are then optimized to reproduce the be- 
havior of the ab-initio bands near the Fermi level. The 
result and the parameters (apart from an overall constant 
energy) of the optimal fit are shown in Fig. |^. 

From the relative weight between the V(l) and the 
V(2) contributions near the Fermi level, p(l)/p(2), we 
learn that there must be a substantial on-site energy 
±eo for the V(l/2) orbitals respectively. We find e = 
0.15 eV. The rung-hopping matrix element t a may be 
expected, on the other hancL,to be quite close to the 
one obtained for a'-NaV^OsEil since the large bending 
angle V(l)-0(1)-V(2) should not substantially affect the 
7r-bonding via the 0(l)-p y orbital. Indeed, we find 
t a = 0.35 eV. We can check whether our estimates for 
eo and t a lead to the correct V(l)-V(2) charge ordering. 
By diagonalizing a simple two-site rung model we find 
the relation 

£0 = P(2)- P (l) (1) 
t a ' ' 2^pOW) ' 

which yields p(l)/p(2) = 2.3 for e /t a = 0.15/0.35, i.e 
p{l) w 0.7 and p(2) « 0.3, in agreement with the DFT 
results. 

We note next, that the crossing along T-Y and Z-U re- 
spectively can be described (within a V-model) only by 

a considerable t% , as_Lt has been noted previously^. 
Recently Yaresko et alxB have discussed that this ma- 
trix element arises naturally when one integrates out the 
coupling between two leg-oxygens of two adjacent edge- 
sharing VO5 pyramids. Our results = — 0.05eV and 
H = -|O n 02eV are close to the values obtained for a'- 
NaV 2 O s B. 

A substantial diagonal hopping matrix element t^ is 
needed in order to explain the fact that the disper- 
sion along T-Y of the lower four V-bands has an opna- 
site behavior with respect to the upper four V-bandsEa, 
as has been pointed out independently for a'-NaX^Os 
by Yaresko et al&B. The substantial contribution of 
td = 0.10 eV pan-.be explained by effective V-O-O-V ex- 
change pathala'cJ. The hopping matrix element along 

the V(l) leg, = — 0.06eV has a similar value to the 

one for a'-NaV205, though the sign of the effective V(2)- 

(2) 

leg hopping parameter, t^ = 0.02 eV, is opposite to the 
expected one. 

(2) 

The result for t\ = 0.05 eV might have been expected 
since its corresponding value for a'-NaX^Os is smallnil. 
The result t± = — 0.18eV is, on the other hand, sub- 
stantially larger and needs some explanation. In terms 
of the bandstructure, tj is determinedpDredominantly 
by the large splitting of the V-d bandso at X and n T, 
(see Figs. || and |J). Also, it has been noted previously^ 
in the context of a'-NaV205 that the bare two-center 
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Slater-Koster matrix elements contributing to t\ can be 
as large as — 0.3eV. The effective t\ is reduced from 
the bare two-center matrix element by the rotation of 
the V-3d orbitals about the crystallogrjaphic y-axis and 
by interference from three-center termsu£3. The contri- 
bution from V-O-V exchange paths depends strongly on 
the relative positions of the three atoms. We have per- 
formed for 7-LiV2C>5 a Slater-Koster analysis and found 
that due to interference effects in between the respective 
dda and ddm contributions, the V(2)-V(2) matrix element 

(2) 

contributing to t\ is smaller than the V(l)-V(l) matrix 
element contributing to t± . In addition, the exchange 
along V(l)-0(4)-V(l) contributing to was found to 
be substantially larger than the exchange V(2)-0(5)-V(2) 

(2) 

contributing to t\ . The particular V-O-V distances and 

angles lead therefore to different t^-p and t±. 

Microscopic model- The microscopic model corre- 
sponding to the results shown in Fig. |5| is that of a spin-^ 
antiferromagnetic Heisenberg chain where the magnetic 
moments are associated dominantly with the sites of the 
V(l) ions in 7-L1V2O5, the contribution is p(l) m 0.7. 
The contribution of about p{2) f» 0.3 of the V(2) to the 
magnetic moment on a V(2) — 0(1) — V(l) rung has 
nevertheless important consequences for the underlying 
microscopic model. For negligible values of p(2) the mi- 
croscopic model could be considered as that of a zig-zag 



chain with a J ~ 4 



since t 



(i) 



-0.06 eV is smaller 



than t 



(i) 



U "b 
0.18eV. In the presence of a non-negligible 



value of p(2) the effective hopping matrix element i[ e ^' 
in between two asymmetric rung states along b is 



4 eff) =p(l)4 1) +p(2)4 2) -2 v ^(lM2)t d 



-0.127eV 



suggesting an asymmetric ladder model. Then, using the 



expression Jb = 2 
transfer gap is E, 



- , and assuming that the charge- 
s' 0.7 eV as in a'-NaVaOjl the ex- 
change integral is Jb w 540 K which overestimates the 
experimental value J exp = 308 K. The degree of charge 
ordering has therefore a substantial influence on the na- 



ture of the magnetic couplingsc-3. 

Conclusions.- We have presented an analysis of DFT 
band-structure calculations for 7-L1V2O5. We find that 
the degree of charge ordering has a substantial influence 
on the nature of the magnetic state. Our results indi- 
cate incomplete charge ordering and 7-LiV20s could in 
this case be viewed as a spin-i asymmetric quarter-filled 
ladder compound. This model would explain the spin 
wave excitation spectri-mi obtained by inelastic neutron 
scattering experimentsEa. Finally, we observe that small 
distortions in the lattice may have substantial effects on 
the interladder V-V hopping matrix element t%. 
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FIG. 1. Crystal structure of (a) 7-L1V2O5 and (b) 
a'-NaV^Os projected in (xz) plane. For 7-L1V2O5 the (xz) 
cut through one of the two equivalent xy planes is shown here. 
The large circles are the V-ions, black and white for V(l) and 
V(2) respectively in 7-L1V2O5 and black for a'-NaV^Os. The 
oxygens are represented by the smaller circles. The alkali-ions 
(Li,Na), shown by grey circles, are located in between the 
planes, close to the bridging oxygens. 
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FIG. 2. Hopping parameters used for the tight-binding 
model for 7-L1V2O5. The arrangement of the V-ions in an 
xy-plane is topologically identical to the one in a'-NaV2 0s, 
i.e. they form a Trellis lattice. The td and the tj are shown 
only partially. Not shown are the onsite energies ±£0 of the 
1/(1/2) sites. 
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FIG. 3. LDA results for LIV2O5. The path is along 
r=(0,0,0), Y=(0 )7 r,0), S=(n,n,0), X=(n,0,0), F, Z=(0,0,7r), 
U=(0,7r,7r), R=(7r,7r,7r), T=(7r,0,7r), Z. The V(l)-3d xy charac- 
ter of the bands is shown with bigger circles. 




FIG. 4. Electron density of the V-3d bands seen in a (xz) 
cut. Shown is the isoline 0.05 e/A 3 . The bigger lobes cor- 
respond to V(l)-3d orbitals and the smaller ones to V(2)-3d 
orbitals. Note that the tilted orbitals point towards the bridg- 
ing oxygens (compare with Fig. |l|) . 
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FIG. 5. Comparison of the tight-binding fit (solid lines) 
with the DFT bands (filled circles). The parameters used 
(see Fig. |) are (in eV) e = 0.15, t a = 0.35, t d = 0.10, 
t[ 1} = ~0.06, 4 2) = 0.02, t[ r> = -0.18, 4 13 = -0- 05 > 
4 2) = 0.05, 4 2) = -0.02. 
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